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Abstract 

Direct computation of the transition time between neighbor resonances in 
the standard map, as a function of the perturbation parameter K, allows for 
improving the accuracy of the critical perturbation value up to K cr — K g < 
2.5 x 10 -4 that is by a factor of about 50 as compared to the previous result 
due to MacKay and Percival. 



As is well known by now a typical structure of the phase space of a few-freedom 
nonlinear dynamical system is characterized by a very complicated admixture of 
both chaotic as well as regular (integrable) components of motion (the so-called 
divided phase space, see, e.g., JI|, |], ||). Statistical properties of such a motion 
are very intricate and unusual. One of the most interesting (and important for 
many applications) problem is the conditions for transition from a local (restricted 
to relatively small regions in phase space) to the global chaos covering the whole 
available phase space. The most studied model of such a transition is described by 
the so-called (canonical) standard map (for history of this model see Q): 

K 

y = y — — ■ sm(2nx) , x—x+y (1) 

where K is the perturbation parameter. In this simple model the transition to global 
chaos corresponds to some exact critical value K = K cr . For K > K cr the motion 
becomes infinite (in momentum y) for some initial conditions while for K < K cr all 
the trajectories are confined within a period of map (|l|): Ay = 1. 

The first idea how to solve this difficult problem was due to Greene |J. First, 
he was able to solve a much simpler problem of the critical perturbation K(r) at 
which a particular invariant Kolmogorov - Arnold - Moser (KAM) curve with the 
rotation number r is destroyed. Critical function K(r) is extremely singular with 
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big dips at everywhere dense set of rational r values (see, e.g., ||). The physical 
mechanism of this behavior (known since Poincare) is explained by resonances in the 
system ([I]) as the rotation number is the ratio of oscillation/perturbation frequencies. 
Whence, the main Greene's idea: to find the 'most irrational' r = r g which would 
correspond to the motion 'most far-off' all the resonances. The former is well known 
in the number theory: r g = [111...] = (\^5 — l)/2 where the first representation is 
a continued fraction. This 'golden' curve was found to be critical at the parameter 
K = K g = 0.97163540631... p, [7|. It was conjectured that for K > K g all invariant 
curves are destroyed @], that is K cr = K g . 

The 'most -irrational' assumption - as plausible as it is - remains a hypothesis. 
The main difficulty is here in that the resonance interaction and overlap, destroy- 
ing invariant curves, depend not only on the resonance spacings, which are indeed 
maximal for r = r g , but also on the amplitudes of those which are not simply 
an arithmetical property. Another argument, based on the analysis of the critical 
function K{r) || ||, also does not prove this principal hypothesis. 

A different approach to the problem - the so-called converse KAM theory - was 
developed in JL0| [TTj] . It relies upon a rigorous criterion for the absence of any 
invariant curve in a certain region. Unfortunately, this criterion can only be checked 
numerically, and besides it provides the upper bound K+ r only (the lower bound 
K~ = K g ). The remaining gap, or the accuracy of K cr : 

(AK) cr = K+ - K cr (2) 

can be made arbitrarily small at the expense of computation time tc which scales 



as 
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t c oc {K+ - K cr y p (3) 



Facing this difficulty, it is natural to recall the first method for calculating the 
critical perturbation used in [|IJ. The method was based on the direct computation 
of trajectories for different K —>■ K g . The criterion of supercriticality of a particular 
K value was very simple: the transition if only a single trajectory in one of two 
neighbor integer resonances (y r = mod 1) through the destroyed critical curve. 
With the computers available at that time the minimal K — 1 has been reached 
only which corresponds to the uncertainty (AK) min = K min — K g = 0.0284. This 



may be compared to the later result (AK) min = 0.0127 (T^ . 

Remarkably, the dependence of the average transition time on parameter K was 
found to be similar to scaling (|3]): 

Fitting three unknown parameters gave: A = 103, p = 2.55, and K cr = 0.989. The 
latter result was rather different from the present value K cr w K g , again because of 
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AK = K-K, 



Figure 1: Direct computation of K cr in standard map: average transition time 
through the destroyed critical curve vs. supercriticality. Circles show numerical 
results for N tr = 400; stars represent 3 single-trajectory runs, including one with 
the minimal AK (^); straight line is relation ([|) with parameters @ fitted from 15 
left-most points (circles). 
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the computation restrictions mentioned above: K > 1, t < 10 7 iterations. Never- 
theless, the fitting Eq.(^) provided a less uncertainty (AK)f = Kf — K g = 0.0174 
as compared to the result from the minimal K. The same is true for data from [|Kj 
where (AK)f = Kf — K g = 0.00236. The latter value was apparently obtained by 
the direct fitting the relation @. Fitting in log-log scale provides a much better 
result: (AK) f = K f - K g = -0.000128 ± 0.000288 that is the remainig uncertainty 
reduces down to 0.000288. 

In both cases the fitted value for the critical perturbation K cr is only true up to a 
certain confidence probability while the minimal K is an exact result: K+ r = K min . 

In the present paper the studies are continued with much better computers. 
The main result is farther considerable increasing of the accuracy (AK) cr . 

To reduce the computation expenses, the transition time was calculated for a 
number of trajectories N tr started near the unstable fixed point of a half-integer 
resonance (y r = 1/2 mod 1), and then run until each of them crossed over to a 
neighbor integer resonance. 

The minimal K value is determined already by the first trajectory escaped from 
the half-integer resonance. In this way the minimal uncertainty 

(AK) min = K min -K g = 0.00025 (5) 

has been achieved with the escape time £ ^ 6.77 x 10 11 itterations which took about 
72 hours of CPU time on ALPHA-4100 computer (see Fig.l). 

The average transition time was computed from N tr = 400 trajectories for each 
of 100 values of K in the interval: 0.0035 < K — K g < 0.35. This costed 36 hours 
of computation. The results are shown in Fig.l. In the whole interval of AK the 
dependence < t(K) > is not exactly a power-law. It becomes so asymptotically for 



K — > K cr as expected from the theory [12]. For this reason, only few smallest K 



values of the function < t(K) > were taken for the final fitting which is also shown 
in Fig.l by the solid line. It is obtained from the fitting 15 left -most points (just 
up to the first big fluctuation) in log-log scale, and corresponds to the following 
parameters in Eq.(|4]): 

(AK)f = 0.000125 ±0.000267, p = 2.959 ± 0.0771, A = 33 ± 8 (6) 

The fitting relative accuracy rms = 0.071 is close to, but somewhat larger than, the 
standard rms = N tr = 0.05. This is seen from the data of 3 single trajectories 
in Fig.l, too. Notice also 2 very big deviations for the average over 400 trajectories 
which nature remains unclear. Interestingly, the relative fitting accuracy of the data 
is considerably higher: rms = 0.02. This would require as many as about 5000 
trajectories in the present method. However, it does not mean that the computation 



of the procedure in |T0[ would be shorter. 

The most important parameter in (|6D is ( AK) f which is zero within statistical 
errors. This farther confirms the Greene hypothesis K cr = K g . The exponent p is 
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also equal to the theoretical value pth = 3.011722 [12| to the fitting accuracy. The 
present value of parameter A is much less than in because of a different (shorter) 
transition between resonances used. The summary of all results is presented in the 
Table below. 



Table. Accuracy of K cr in standard map 





(AK) flt 


Reference 


exact 


probable 








2.84 x 10" 2 


1.74 x 1(T 2 


B 


1.27 x 1CT 2 


3.36 x ICT 3 
±1. x ICT 3 




H3 






-1.28 x ICT 4 




iq 






±2.88 x ICT 4 


our fit 


2.5 x ICT 4 


1.25 x ICT 4 


present 




±2.67 x 1(T 4 


paper 



A serious difficulty in such a numerical approach to the problem is the computa- 
tion accuracy. This was mentioned also in |H| but no estimate for the computation 
errors was given, apparently because of a very complicated numerical procedure. 
Even in a much simpler method [Q, accepted in the present study, the effect of 
noise turned out to be rather complicated. Special numerical experiments were 
done to clarify the question. To this end, a random perturbation of amplitude v 
was introduced in both equations (|l]). The results are shown in Fig. 2. 

Typically, the transition time becomes less than that without noise, and saturates 
below some critical noise-dependent value of K: AK ^ B(u). However, in some 
cases the average transition time considerably grows, as an example in Fig. 2 demon- 
strates, apparently due to a sharp increase of the fluctuations near the crossover 
from normal (noisefree) dependence of < t(AK) > to the saturation. In turn, these 
fluctuations are apparently explained by the noise-induced diffusion into some of 
many small domains of regular motion within the critical structure. 

A rough estimate for unknown function B(u) can be obtained as follows. The 
transition time is primarily determined by the width 5y ~ (AK) 2 of the chaotic 
layer around destroyed critical curve [12|, OL |3|] while the diffusion time through this 
layer t ~ l/AK |14|, [13|, |15|, |I6| . Noise decreases this time down to t u ~ (5y) 2 /u 2 . 
Hence, the crossover corresponds to t u ~ t , whence: 

B(v) w A ■ v h (7) 
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Figure 2: The effect of noise on the supercritical transition time: straight line is 
the fit of the noisefree computation results (cf. Fig.l); points connected by lines 
represent the impact of noise with amplitude v computed for N tr = 10; numbers at 
lines are -log(z/) values (logarithm decimal). 
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Figure 3: Noise scaling: circles give the empirical crossover values B(v) as a function 
of noise amplitude v connected and extrapolated by the straight line (0); the upper 
horizontal dotted line shows the minimal AK in computation with noise while the 
lower line indicates (AK) min (H) achieved in the main double-precision computation 
(see Fig.l) with the accuracy roughly corresponding to log(V) ~ —15; all logarithms 
are decimal. 
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with b = 2/5. Fitting the empirical data in Fig. 2 in log-log scale gives: b pa 0.39 ± 
0.012, which is surprisingly close to the theoretical estimate, and A pa 0.9716 ±0.054 
(Fig. 3). The fitting accuracy is also fairly good: the relative rms = 0.019. Moreover, 
below crossover (AK < B{y)) the width Sy as well as the diffusion time depend on 
v only, and hence the transition time remains approximately constant for a given 
v (Fig. 2). In any event, the minimal (AK) m i n (|5|), which is the main result of 
the present study, is well above the expected limitation for the double-precision 
computation (see Fig. 3). 

In conclusion, the direct approach a la to the problem of the critical pertur- 
bation in the standard map does further confirm Greene's hypothesis K cr = K g with 
a much better exact upper bound ©: K cr — K g < 2.5 x 10~ 4 . 

Still another recent confirmation of this conjecture (curiously, with roughly the 
same statistical accuracy @) has been inferred from a detailed study of the critical 
structure at the chaos-chaos border in standard map for K = K g fTE|| . 
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